BAR4R-CONF-04/20 
SLAC-PUB-10610 
August 2004 



Measurement of the Branching Fractions for the Decays 
^0 _^ ^*-^_^+^ ^0 _^ ^0 _^ D*^pp and D^pp 



The BAR4R Collaboration 
February 7, 2008 



Abstract 

The meson decay modes D*~ppiT^ , D^pp-K^ , D*^pp, and B^ D^pp are studied in 

a sample of 124 x 10^ BB pairs collected with the BABAR detector at the PEP-II colhder. The 
decay of B'^ D^ppn^ is observed for the first time, with a measured branching fraction B{B^ 
D~ppTT^) = (3.80 lb 0.35 lb 0.46) x 10~^. The following branching fractions are also determined: 
B{B^ D*-ppn+) = (5.61 ± 0.59 ± 0.73) x 10"^, B{B^ D*^pp) = (0.67 ± 0.21 ± 0.09) x 10"^, 
and B{B^ D pp) = (1.24ib0.14ib0.12) x 10"'^. In each decay mode the invariant mass spectra of 
the charmed mesons and the baryons are compared with a pure phase-space hypothesis in order to 
gain insight into the B meson decay dynamics. In particular, the Dalitz plots of D*^p versus D*^p 
for B^ — > D*^pp and of D^p versus D^p for B^ D^pp are presented. All results are preliminary. 
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1 INTRODUCTION 



The very successful performance of the two B factories PEP II and KEKB enables the study of 
B meson decays with unprecedented sensitivity. In this paper B meson decays to final states 
which include a charmed meson and a baryon anti-baryon pair are studied. The observation of 
the B^ — > D*~ppTr~^ and B° D*- 'pn decays by CLEO and the color-suppressed B^ ^ D pp 
and B^ D*^pp decay modes by Belle suggest the dominance of multi-body final states in 
decays of B mesons into baryons 3 . In this paper we present the measurements of the branching 
fractions for the following four decay modes: B^ D*~ppK~^ , B° ^D- B^ ^ D*' pp, and 

pp. Observation of modes proposed here would help clarify the dynamics of weak decays 
of B involving baryons jlj . 

Since the branching fractions of multi-body decays are sizable [3], it is natural to ask whether 
such final states are actually the products of an intermediate two-body channel. If this is the case, 
then these initial two-body decays may involve a baryon-antibaryon bound states (NN) [HI Ej or a 
charmed pentaquark {DN, which can be cuudd for example) [SI E] or a (nonexotic) heavy charmed 
baryon. Motivated by these considerations, the invariant mass spectrum of the baryon-antibaryon 
and the invariant mass spectra of the charmed meson and baryon are investigated. In particular, 
the Dalitz plots for D*^p versus D*^p and D^p versus D^p for the B^ — > D^pp and B^ D*^pp 
decay modes, respectively, are presented. The inclusion of the charge conjugate modes is implicit 
throughout this report. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage 
ring. The sample corresponds to an integrated luminosity of 112.5 fb~^ at the T(45') resonance 
(on-resonance) and 11.9 fb~^ taken 40MeV below the T(4S') resonance (off-resonance). The on- 
resonance sample contains about 124.1 x 10^ BB pairs. The collider is operated with asymmetric 
beam energies, producing a boost of (3"f ~ 0.56 of the T{AS) along the collision axis. 

The BABAR detector is optimized for asymmetric energy collisions at a center-of-mass (CM) 
energy corresponding to the T{AS) resonance. The detector is described in detail in reference |1()| . 
Charged particle tracking is provided by a five-layer double-sided silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) contained within the magnetic field of a 1.5T superconducting 
solenoid. The tracking system provides momentum reconstruction of charged particles and measures 
energy-loss (dE/dx) for particle identification. Additional charged K-tt particle identification is 
provided by a ring- imaging Cherenkov detector (DIRC), which exploits the total internal refiection 
of Cherenkov photons within synthetic quartz bars. The energies of neutral particles are measured 
by an electromagnetic calorimeter (EMC) composed of 6580 CsI(Tl) crystals. The magnetic flux 
return of the solenoid (IFR) is instrumented with resistive plate chambers in order to provide muon 
and neutral hadron identification. 

A GEANT4-based JAj Monte Carlo (MC) simulation is used to model the signal efficiency 
and the physics backgrounds. Simulation samples equivalent to approximately three times the 
accumulated data were used to model BB events, and samples equivalent to approximately one 
times the accumulated data were used to model e+e" uu, dd, ss, and cc events. 
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3 ANALYSIS METHOD 



The meson is reconstructed in the following four decay modes: i? D* ppTT~^, 
B° ^ D pp, and B^ D^pp. The D*' and D*^ are reconstructed as D tt and D^n^, respec- 
tively. The D~ candidates are reconstructed in the decay mode K~^tt~tt~. The candidates are 
reconstructed in the K~^7r~ mode. 

3.1 SELECTION 

The B^ reconstruction proceeds as follows. First we reconstruct the and candidates in the 
decay modes noted above. For all decay modes, charged kaons are distinguished from the pions and 
protons with energy- loss (dE/dx) information in the tracking system and the Cherenkov angle and 
the number of photons measured by the DIRC. For the B^ D*^pp decay mode only, the pion 
from K^7r~ decay must not be identified as either an electron, muon, kaon or proton. The 

mass of the reconstructed or D~ candidates must be within 3 standard deviations of the fitted 
mean of reconstructed or D~ mass for each decay mode. Finally, the daughter tracks from the 

or D~ are required to be consistent with originating from a common vertex. 

For the decay modes B^ D*^ppT:^ and B^ — > D*^pp, D*~ and D*^ candidates are recon- 
structed by combining the candidate with a soft track or a soft n^, respectively. The momentum 
of the soft track or soft vr^ in the CM frame must be less than 0.45 GeV/c. The ir^ candidate is re- 
quired to have a two-photon invariant mass between 0.116-0.150 GeV/c^, and its daughter photon 
candidates must have a minimum cluster energy of 30 MeV. The mass difference between D* and 

(AM = mo* - rnj^o) is required to be less than 0.1486 GeV/c^ and 0.1453 GeV/c^ for D*' 
and D*^ candidates, respectively. The AM requirement removes crossfeed from B^ — > D*^pp into 
B^ D^pp. For each track except for the soft pion from D*~ decays, the transverse momentum 
Pt must be larger than 0.1 GeV/c in order to improve the quality of the vertex fit. 

To reconstruct the candidate B^ D*^pp-K^ and B^ D^ppK~^ decay modes, both proton 
and anti-proton candidates are distinguished from pions and kaons on the basis of energy-loss 
{dE/dx) information in the tracking system and the Cherenkov angle and the number of photons 
measured by DIRC. The reconstructed D*~ and D~ , respectively, are combined with an identified 
proton and anti-proton pair and a track. We require that the pion candidate track must not be 
identified as either an electron, a muon, a kaon or a proton. For the decay modes B^ D* pp and 
pp, the D*^ and candidates, respectively, are combined with an identified proton and 
anti-proton pair. All the daughters of the B^ must originate from a common vertex. 

Two additional kinematic variables are used to identify the reconstructed B^ candidates |10j . 
The first is the beam-energy-substituted mass, rriES = [{^'cm/'^ + Pi ' Pb)'^/^? ~ p'b]^^'^^ where 
EcM is the total CM energy of the e~^e~ collision, (Ei, pj is the four-momentum of the initial e~^e~ 
system and is the momentum of the reconstructed B candidate, both measured in the laboratory 
frame. The second variable is AE = E'^ — Ecm/'^i where E"^ is the i?-candidate energy in the CM 
frame. The B meson candidates are defined by requiring: 5.20GeV/c2 < mss < 5.29 GeV/c^ and 
lAi^l < 0.12 GeV. These two variables are then used in a maximum likelihood fit to extract signal 
and background yields. 

If more than one suitable B^ candidate is reconstructed in an event, then for each B^ decay 
mode one best candidate is selected. The best candidate selection algorithm is based on a x^-like 
quantity constructed from the difference between the D mass and/or mass difference AM for the 
candidate and the nominal value |12J. For a given event, the candidate with the lowest value of 
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is selected for each decay mode. For the D* ppir^ and B^ — > D* pp decays, the best 

candidate is selected based on the mass of the candidate and AM. For the B^ D~ppTT~^ and 
pp decay modes, the reconstructed D meson mass is used for best candidate selection. 
To suppress background from two-jet-like e~^e~ qq continuum, variables that characterize 
the event topology are used. We require cos Ot^r < 0-9, where Othr is the angle between the thrust 
axis of the B^ candidate and that of the rest of the event. This requirement eliminates 63% of the 
continuum background and retains 88% of the signal events. For further continuum background 
suppression we require that the ratio of the second to zeroth Fox- Wolfram moments ^H] is less than 
0.35. 

The signal efficiency in each B decay mode after applying all selection criteria is shown in 
Table ^ The efficiencies listed in the table are obtained from signal MC simulation assuming pure 
phase-space for the B decay model. 

Table 1: Summary of selection efficiencies obtained from MC simulation based on a phase-space 
model and fitted yields from data for the B decay modes in this report. The error on the efficiency 
comes from the Monte Carlo statistics. The statistical significance is calculated as -v/2Alog£, where 
A log C is the log-likelihood difference between a signal hypothesis corresponding to the yield and 
that corresponding to a null yield. 



Mode 




Efficiency, % 


Signal yield 


Statistical significance 


B^^ 


D* ppir^ 


6.97 ± 0.10 


130±14 


18CT 


B^^ 


D^ppTT^ 


5.87 ± 0.10 


238±22 


17cj 


B^^ 


D pp 


5.53 ± 0.09 


13± 4 


So- 


B^^ 


D pp 


16.53± 0.15 


96 ± 11 


17CT 



3.2 THE MAXIMUM LIKELIHOOD FIT 

We perform an unbinned extended maximum likelihood (ML) fit to extract the signal and back- 
ground yields. The variables mES) are used to separate signal from backgrounds. The data 
sample is assumed to consist of signal and combinatorial backgrounds, which arise from random 
combinations of charged or neutral candidates from both continuum and BB events. The extended 
likelihood for a sample of N candidates is 

N 

C = e-^' ■ X{{N,,g ■ Vs^g{mESi, AEi) + N^tg ■ VbkgimESi, A^OI (1) 

i=l 

where Nsig and N^^g are the number of signal and background yields (to be determined by the fit), 
respectively, and A^' = Nsig+Nf^kg- The probability density functions (PDFs) Vsig and Vbkg are 
the product of the PDFs of two discriminating variables. The signal PDF is thus given by Vgig = 
V{raESi) ■ 'Pi^Ei). The signal PDFs are decomposed into two parts with distinct distributions: 

• Signal class I: signal events that are correctly reconstructed or signal events that are misre- 
constructed due to a random slow vr" or vr*^ assigned to a D*^ or D*^ decay. Since the rriEs 
and AE distributions for these two kinds of events do not significantly differ from each other, 
we put them into one signal class in the ML fit. 
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• Signal class II: signal events that are misreconstructed due to wrong tracks from B direct 
decays or from a wrongly reconstructed D meson. The ttt-es and AE distributions are sig- 
nificantly different from those of signal class I. Hereafter signal class II events are also called 
Self-Cross-Feed (SCF) events. The fraction of this kind of misreconstructed signal event is 
estimated by MC simulation. 

The PDFs of mss and AE for signal events are taken from the MC simulation, with the exception 
that the means of the signal Gaussian for the tit-es and AE PDFs are free to vary in the fit. The 
niES and AE PDFs for the combinatorial backgrounds, which include B background and continuum 
background, are described by two free parameters. One of these is from an ARGUS function |14j 
used to describe the tties shape, the other is from a first-order polynomial for the AE shape. A 
total of six parameters, including signal and background yields and all the parameters related to 
background PDFs, are varied in the fit. 



4 PHYSICS RESULTS 

The signal yield from the ML fit and statistical significance for each of the B modes are summarized 
in Tabled The statistical significance is calculated as ^/2ATogC, where A log C is the log-likelihood 
difference between a signal hypothesis corresponding to the yield and that corresponding to a null 
yield. Figure H shows projection plots of rriEs and AE from the likelihood fit. The fit shows 
satisfactory agreement with distributions of the discriminating variables. 

To correct for the efficiency variation across the Dalitz plot, which is of order 20% for B^ — > 
D*^pp and 5% for B^ — > D^pp three-body decays, each event is assigned a signal weight, Wgig, 
which takes correlations into account: 

where A'^- is the number of events, V{m^sT is the PDF of the j*^ component given the event 
values of mEs and AE, and Vsig,] is the signal row of the covariance matrix of the components 
yields obtained from the likelihood fit. We determine the efficiency as a function of position on the 
Dalitz plot from simulated signal events. The branching fractions of B^ D*^pp and B^ — > D^pp 
three-body decays are corrected by the Dalitz plot dependent efficiencies. The effect of the D*~ 
polarization in the B^ — > D*^pp mode is neglected in this Dalitz plot dependent efficiency correction. 

We estimate the efficiencies for 4-body decays with the same method as for three-body decays. 
We find that the differences are 7.5% and 3.1% from the efficiencies obtained from a pure phase- 
space MC for the B^ — > D*~pp7r~^ and B° D' "ppTT^, respectively. We will assign 7.5% and 3.1% 
as additional systematic errors for the two four-body decays, respectively. 

In each of the four B decay modes, the signal region is defined as: —0.035 < AE < 0.03 GeV and 
w-ES > 5.27 GeV. The ttt-es distributions after applying the requirement —0.035 < AE < 0.03 GeV, 
and AE distributions after applying the requirement rriES > 5.27 GeV/c^, are shown in Fig. IJJfor 
all four decay modes. We also compare the invariant mass spectra for charmed meson and baryon 
combinations with a pure phase space hypothesis for B decays. Figure 01 shows the invariant mass 
of D*~p and D*~p for signal candidates in the B^ — > D*~ppTT~^ mode. The open histogram is the 
expected distribution from B^ D*^ppiT~^ phase space signal MC simulation which is normalized 
to the fitted signal yield in data. The cross-hatched histograms describe the contributions from 
different background processes, which is normalized to the fitted background yield in data. 
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m,, (GeV) 



AE (GeV) 



Figure 1: mss (left-hand column) and AE (right-hand column) distributions in data. The plots 
from top to bottom are for D*^pp7r^, D^pp-K^, 5° D*^pp, and 5° Wpp, 

respectively. The solid curve represents a projection of the maximum likelihood fit. The dashed 
curve represents the contribution from combinatorial background. 
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Figure 2: mss (left-hand column) distributions from data after applying requirement —0.035 < 
AE < 0.03 GeV. The AE (right-hand column) distributions from data after applying requirement 
TTT-ES > 5.27 GeV/c^. The plots from top to bottom are for S° — £)*~pp7r+, D~pp7r+, 
pp, and 5° — > D^pp, respectively. 
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M(D p/D '*p) (GeV) M(d' p/D '+p) (GeV) 



Figure 3: The invariant mass distribution of D*~p (left) and D*^p (right) combinations from 
signal candidates in the — > D*~'pp-K^ decay mode. The cross-hatched histograms describe the 
contributions from different background processes. The open histogram is the expected contribution 
from 5° D*- 'ppiT'^ phase space signal MC simulation. 




M(D"p/D*+p) (GeV) M(D"p/D**p) (GeV) 

Figure 4: The invariant mass distribution of D~p (left) and D~p (right) combinations from signal 
candidates in the — > D~ppTT'^ decay mode. The cross-hatched histograms describe the contribu- 
tions from different background processes. The open histogram is the expected contribution from 
the B^ — > D~ppTT~^ phase space signal MC simulation. 



The shape of the invariant mass spectra in Fig. El show that the D*~p distribution in data is 
only marginally consistent with a phase space model, while the D*~p is quite consistent with phase 
space. A similar observation holds for the D^p and D^p distribution from signal candidates in the 
mode B^ D^pp-n^ as shown in Fig. 

As shown in Fig. El the shapes of the invariant mass spectra from a phase-space model for the 
D*^p and D*^p distributions in the B^ D*^pp mode, are consistent with data given the current 
low statistics. However, the shape for the D^p distribution in the B^ D^pp mode is inconsistent 
and it appears to deviate from the phase-space expectation, as illustrated in the left-hand plot 
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Figure 5: The invariant mass distribution of D*^p (left) and D*^p (right) combinations from signal 
candidates in the — > D*^pp decay mode. The cross-hatched histograms describe the contribu- 
tions from different background processes. The open histogram is the expected contribution from 
the — > D*^pp phase space signal MC simulation. 




M(i5"p/D"p) (GeV) M(d"p/D"p) (GeV) 



Figure 6: The invariant mass distribution of D^p (left) and D^p (right) combinations from signal 
candidates in the B^ D^pp decay mode. The cross-hatched histograms describe the contributions 
from different background processes. The open histogram is the expected contribution from the 
pp phase space signal MC simulation. 



of Fig. IHI For D^p, the distribution in data is consistent with the phase space, as shown in the 
right-hand plot of Fig. |HJ 

To further compare the observed mass distributions in data with uniform phase space, we 
show the Dalitz plot distributions for the two three-body decay modes in Figures [7| and |S1 The 
distribution of m'^{D*^p) versus m'^{D*^p) for signal candidates in the B^ D*^pp mode is shown 
in Figure [3 the distribution between m?{D^p) and m?{D^p) for the B^ D^pp mode is shown in 
Figure IHl The triangles and dots correspond to data and the phase space MC simulation. Figure [3 
shows a Dalitz plot for the mode B^ D*^pp. The comparison with the phase space distribution 
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8 10 12 14 16 18 

M^(D*°p/D*V)(GeV^) 

Figure 7: Dalitz plot for the mode D*^pp. The horizontal and vertical axes correspond to the 

squared invariant masses of the D*^p and D*^p, respectively. The triangles and dots correspond to 
signal events in data and the phase space simulation, respectively. 




8 10 12 14 16 18 



M\D p/D^pXCeV^) 

Figure 8: Dalitz plot for the mode — > D^pp. The horizontal and vertical axes correspond to the 
squared invariant masses of the D^p and D^p, respectively. The triangles and dots correspond to 
signal events in data and the phase space simulation, respectively. 
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Figure 9: Invariant mass distributions for pp pairs in the signal sample for (a) for D*~pp7r~^] 
(b) for D~ppTT~^\ (c) for B^ — > D*^pp; and (d) for B^ — > D^pp, respectively. The cross- 
hatched histograms describe the contributions from different background processes. The open 
histogram is the expected contribution from the phase space signal MC simulation. 
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is limited by low statistics. The equivalent plot for the mode — > D^pp, shown in Figure |HJ shows 
a deviation from the pure phase space model. In addition, the Dalitz plots display a threshold 
enhancement in the pp invariant mass spectrum. Figure shows the invariant masses of pp for each 
B decay mode. A threshold enhancement is observed for each B decay mode. A similar effect has 
also been observed by other experiments I16j. 

The compatibility of data and MC simulation was quantified by means of a Kolmogorov-Smirnov 
test. Table m shows the probabilities of compatibility of data and MC simulation, where the shape 
of the data points and the MC histograms shown in Figures f.SI6l andl9l were compared. The first and 
second columns correspond to the invariant mass spectra of the relevant D-meson and p combination 
(the corresponding minimal quark content is given in parentheses). Charge conjugation is assumed. 
The last column corresponds to the invariant mass spectra of the pp system. We draw a number 
of observations from the calculated probabilities in Table |21 For the first column, the top two 
rows are marginally consistent with phase space, while the third is consistent, and the last row is 
inconsistent. The entire second column is consistent with phase space, while in the third column 
only the third row is consistent, and all others are inconsistent. 

Table 2: Probabilities of obtaining the observed result under the null hypothesis that the data is 
sampled from the phase space model. 



B° D*-ppTT+ 


D* p{cduud) 
1.32x10-2 


D* p{cduud) 
0.83 


M{pp) 
1.22x10-^ 


BO D~pp7T+ 


D^p{cduud) 
1.78x10-2 


D^p{cduud) 
1.00 


M{pp) 
2.58x10"'^ 


B° D*^pp 


D*^p{cuuud) 
0.13 


D*^p{cuuud) 
0.68 


M{pp) 
0.32 


B^ D^pp 


D^picuuud) 
3.51x10"'' 


D^p{cuuud) 
0.77 


M{pp) 
3.26x10-^ 



5 SYSTEMATIC STUDIES 

Systematic uncertainties in the ML fit originate from the modeling of the PDFs. We vary the PDF 
parameters within their respective uncertainties, and derive the associated systematic errors. The 
SCF fractions can be floated in the fit to data for the decays B^ — > D*~ppTT~^ and B^ — > D~ppK~^ . 
The differences between MC and data are 8% and 16% for the two decays, respectively. For the 
systematic uncertainty due to the SCF fraction, we vary the SCF fraction by ±10% in the fit 
according to the above validation from data. We also perform fits to large MC samples with signal 
and combinatorial background. No bias is observed in these tests. 

The systematic errors in the efficiencies are for the track finding (1.2% - 1.4% per track and 
2.2% for the slow charged vr), particle identification (0.1 - 0.3 %) and vr'^ reconstruction (5.6%). 
The systematic error in each B decay mode arising from variations of the selection criteria is also 
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Table 3: Summary of the systematic uncertainties for all B modes in this report. 



Error source 




(m 7oj 






Event selection 


Q 1 
O.i 


o.y 


O.D 


1.0 


Signal model 


3.0 


3.2 


9.0 


0.5 


Track finding 


9.2 


8.4 


4.8 


4.8 


vr*^ reconstruction 






5.6 




Particle ID 


0.1 


0.3 


0.1 


0.1 


Phase space MC 


7.5 


3.1 






Number of B mesons 


1.1 


1.1 


1.1 


1.1 


Daughter branching fraction 


2.5 


6.5 


5.3 


2.4 


Total 


13.0 


12.2 


13.9 


9.4 



shown in Table El The reconstruction efficiencies depend on the final state distributions due to B 
decay dynamics. As discussed earlier, we assign 7.5% and 3.1% systematic uncertainties to account 
for the uncertainty on the efficiencies from the phase-space MC simulation for B^ — > D*~pp-K^ and 
B^ ^ D- 'ppiT^ four-body decays, respectively. Other systematic effects are from event-selection 
criteria, daughter branching fractions J2], MC statistics, and the number of B mesons in the 
sample |15| . The contributions to the systematic errors on signal yields are summarized in Table 01 

6 SUMMARY 

The B^ meson decay modes B^ D*~ppTT^ , D^ppn^ , D*^pp, and B^ D^pp have been studied 
in a data sample equivalent 113 fb~^ of integrated luminosity. Table ^ summarizes the branching 
fractions obtained for the decays B^ — > D*~ppTT~^ , B^ ^D- pp and B'^ ^ D pp. 

These results are compared with previous measurements performed by the CLEO ^ and Belle [21 
collaborations, when possible, and good agreement is found. The B^ — > D~pp-K~^ decay has been 
observed for the first time with the branching fraction B{B^ D^ppn^) = (3.80 it 0.35 it 0.46) x 
10~^. All results are preliminary. 

Table 4: The branching ratios (in units of 10"'*) for the B^ modes are measured here. Statistical 
(first) and systematic (second) errors are given. The significance includes both the statistical 
error and the systematic error. Results obtained here are compared with results published by the 
CLEO [U and Belle |2] collaborations (when available). 



Final state 




Significance 


Reference i3(10-'') 


^0 ^ D*-ppTT+ 
B^ D-ppTr+ 
S° ^ D*^pp 
^ D^pp 


5.61 ±0.59 ±0.73 
3.80 ± 0.35 ± 0.46 
0.67 ±0.21 ±0.09 
1.24 ±0.14 ±0.12 


13(7 

13a 
4.O0- 

12a 


i.2oj:[J:i ±0.21 

1.18 ±0.15 ±0.16 
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